Time-delayed nonsequential double ionization with few-cycle laser pulses: importance 

of the carrier-envelope phase 
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We perform theoretical investigations of laser-induced nonsequential double ionization with few 
cycle pulses, with particular emphasis on the dependence of the electron-momentum distributions 
on the carrier-envelope phase. We focus on the recollision-excitation with subsequent tunneling 
ionization (RESI) pathway, in which a released electron, upon return to its parent ion, gives part 
of its kinetic energy to promote a second electron to an excited state. At a subsequent time, the 
second electron is freed through tunneling ionization. We show that the RESI electron-momentum 
distributions vary dramatically with regard to the carrier-envelope phase. By performing a detailed 
analysis of the dynamics of the two active electrons in terms of quantum orbits, we relate the shapes 
and the momentum regions populated by such distributions to the dominant set of orbits along 
which rescattering of the first electron and ionization of the second electron occurs. These orbits 
can be manipulated by varying the carrier-envelope phase. This opens a wide range of possibilities 
for controlling correlated attosecond electron emission by an adequate pulse choice. 
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I. INTRODUCTION 

Ultrashort laser pulses with durations of a few opti- 
cal cycles can reach very high intensities and, still, carry 
much less energy than their longer counterparts 1} . This 
enables a wide range of applications, such as controlling 
the collective electron motion in condensed matter 
manipulating chemical reactions Q, and the generation 
of isolated XU V attosecond pulses 0, @ . In this pulse- 
length regime, the carrier-envelope phase (CEP), i.e., the 
phase of the carrier frequency with respect to the pulse 
envelope, dramatically influences strong-field phenom- 
ena, such as high-harmonic generation (HHG) above- 
threshold ionization (ATI) Q and laser-induced nonse- 
quential double ionization (NSDI) . Specifically for NSDI, 
it has been shown that the electron-momentum distri- 
butions change their shapes considerably when this pa- 
rameter varies. In fact, asymmetric electron momentum 
distributions have been identified both theoretically [7|- 
[l2| and experimentally [l3|-Qj| for NSDI with few-cycle 
driving pulses. Most of these studies have focused on the 
electron-impact ionization pathway, in which the second 
electron has enough energy to overcome the second ion- 
ization potential and both electrons reach the continuum 
simultaneously. 

In recent experiments, however, a similar effect has also 
been observed in a parameter range for which another 
NSDI pathway is prevalent: Recollision excitation with 
subsequent tunneling ionization (RESI) [3, [lj|. RESI 
happens when the first electron, upon its return to the 
core, promotes the second electron to an excited bound 
state, from which, with laser assistance, it subsequently 
tunnels [l6j]. In RESI, the first electron leaves immedi- 
ately upon rescattering, while the second electron tunnel 
ionizes near the subsequent field maximum. Hence, there 



is a time delay between rescattering of the first electron 
and tunnel ionization of the second electron. For a com- 
prehensive discussion of RESI see our recent review [17| . 

The above-mentioned experimental evidence shows 
asymmetric electron momentum distributions, as func- 
tions of the momentum components Pi||,P2|| parallel to 
the laser-field polarization, whose shape varies dramat- 
ically with the CEP. For instance, in [lj| asymmetric 
distributions have been observed, whose probability den- 
sities, depending on this phase, are either stronger in the 
positive half axes p^ = 0, Pji\ > 0, with i.j = 1,2 and 
i =/= j, or in the negative half axes Pi\\ = 0, Pj\\ < 0, 
defined in the parallel-momentum plane Pi\\P2\\- As the 
CEP is varied, the momentum region in which the corre- 
lated probability density is larger shifts from one half axis 
to the other. This behavior resembles to a great extent 
that observed for electron-impact ionization [l3| . which 
has been explained by us in previous work in terms of a 
shift in the dominant set of orbits along which inelastic 
rescattering occurs @, H|. The main difference is that, 
instead of being located either at the positive or nega- 
tive momentum half axes mentioned above, the electron- 
impact NSDI distributions populated either the first or 
the third quadrant of the Pii|P2|| plane. Depending on 
the CEP range chosen, they shifted from the first to the 
third quadrant of the parallel-momentum plane, or vice 
versa. 

In the present work, we address the question whether, 
similarly to what happens in the direct pathway, it is pos- 
sible to associate dominant sets of trajectories to specific 
values of the CEP and to specific shapes of the electron- 
momentum distributions. For that purpose, we apply the 
analytical RESI model developed by us in previous pub- 
lications miHH to few-cycle pulses of different CEPs. In 
this model, the Feynman diagram corresponding to RESI 
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has been considered from the outset, and the pertaining 
transition amplitude has been calculated in the strong- 
field approximation using the steepest descent method. 
A particularly important issue in the context of RESI 
with few-cycle driving pulses is that particular care must 
be taken with regard to causality if the above-mentioned 
methods are used. This issue has been addressed by us 
in a recent publication [27] . 

This article is organized as follows. In Sec. [TT1 we pro- 
vide a brief summary of the method developed in (2lLl2^. 
In Sec. IIII1 we determine the dominant sets of orbits 
for the first and second electron, from the solutions of 
the saddle-point equations. Subsequently, in Sec. IIVI we 
compute the RESI electron-momentum distributions and 
analyze their CEP dependence in terms of such orbits. 
Finally, in Sec. [V] we state the main conclusions to be 
drawn from this work. 



II. MODEL 

A. Transition amplitude 

Within the strong-field approximation, the transition 
amplitude corresponding to RESI reads [21] 
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are the form factors related to the ionization of the first 
electron, recollision of the first electron with excitation 
of the second electron, and tunnel ionization of the sec- 
ond electron [Eqs. ([3]), ((H) and (O, respectively]. These 
form factors provide information about the binding po- 
tential V^(ri) and Vion(r2) "seen" by the first and the 
second electron, respectively, and about the interaction 
Vi2(i"i,r 2 ) of the first electron with the core. 

Eq. ([T]) describes a process in which the first electron, 
initially bound in the ground state l^f 9 ^} with energy 
e[ 9 ^ , is freed by tunnel ionization at time t" into a Volkov 
state \k(t")). Thereafter, it propagates in the continuum 
from the time t" to the time t' with intermediate momen- 
tum k. At a time t' , it rescatters inelastically with the 
core and, through the interaction V\ 2 , excites a second 
electron from the ground state \ip^) °f ^ ne sm gly ionized 
target to the state {ip^)- The energies of the ground and 
excited states of the singly ionized target are E { 2 9) and 

E 2 j respectively. The first electron reaches the detector 
with final momentum pi immediately after rescattering. 
The second electron remains bound until a later time t, 
when it is released by tunnel ionization into a Volkov 
state |p>2 (£))■ It reaches the detector with final momen- 
tum p2. In the above-stated equations, k(r) = k + A(r) 
and p„(r) = p„+A(r) (r = t,t',t") in the length gauge, 
and k(r) = k and p„(r) = p„ in the velocity gauge, with 
n = 1,2. In [21| . we have verified that, in practice, the 
results obtained in both gauges are nearly identical [32j . 



B. Saddle-point equations 

The transition amplitude |T]) is computed using the 
steepest descent method. Hence, one must find the 
values of t,t',t" and k for which the action (J5J) is 
stationary, i.e., for which dt"S(pi, p 2 , k, t, t', t") = 
d t >S(pi,p 2 ,k,i,i',;") = d t S(p 1 ,p 2 ,-k,t,t',t") = 0, and 
d k S(pi, P 2,k,i,i',i") = 0. 

The conditions upon t", k and t' give the saddle-point 
equations 



(ff) 
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and 



[ Pl + A(t') 2 ] = [k + A(t')f - 2(E { 2 9) - E { 2 % (8) 
while the condition upon t leads to 



[ P2 +A(t)] 2 = -2i4 e) . (9) 

Equation © expresses the fact that the energy of the 
first electron is conserved at the instant t" at which tun- 
neling ionization occurs. Equation (J7J) fixes the interme- 
diate momentum of the first electron in order to guaran- 
tee its return to its parent ion. Equation (|5J) states that, 
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upon rescattering, the first electron transferred a fraction 
E cxc = E ( 2 a) - E ( 2 ] of its kinetic energy [k + A(t')f /2 
upon return, and reached the detector with final kinetic 
energy [pi + A(i')r /2. Finally, Eq. © gives the en- 
ergy conservation for the second electron upon tunnel- 
ing, which reaches the detector with final momentum 
P2. Note that Eqs. (j6|) and ((9]) have no real solution, 
since tunneling has no classical counterpart. If written in 
terms of the electron momentum components (p n \\,Pn±)> 
n = 1,2 parallel and perpendicular to the laser- field 
polarization, Eqs. © and © give the kinematic con- 
straints related to the first and second electron, respec- 
tively. These constraints have been discussed in detail in 
Refs. [2l|, |23j and [27j for monochromatic fields and few- 
cycle pulses, respectively, and will be briefly mentioned 
here. 

Explicitly, from Eq. © one obtains the condition 

- A(t') - V / 2AE < pin < -A(t') + y/2AE, (10) 

where AE = E^i n (t', t") — E cxc denotes the energy differ- 
ence between the kinetic energy Ey.i n (t\t") of the first 
electron upon return and the energy E cxc = E cxc +p\ ± /2. 
The latter is an effective excitation energy, which in- 
creases with perpendicular momentum pu_. Inside the 
boundaries defined by Eq. (|10l) , rescattering has a classi- 
cal counterpart and the probability density associated 
with it is significant, while, outside those boundaries, 
rescattering is not classically allowed to occur and the 
corresponding probability density is vanishingly small. 
The largest region will be obtained for vanishing trans- 
verse momentum pu_. 

This implies that (i) the region in the parallel mo- 
mentum plane determined by rescattering of the first 
electron is centered around pin = —A(t'), where A(t') 
is the vector potential at the time the first electron re- 
turns, and (ii) the extension of this region is determined 
by the energy difference AE. In order to draw upper 
bounds for such momentum regions, it is useful to as- 
sume that pij^ = and consider the maximal kinetic en- 
ergy E^ x (t' ,t") the first electron may have upon return. 
If E^ x (t',t") » E cxc , this region will be very large. In 
contrast, if E^^(t', t") ~ E cxc , the energy of the return- 
ing electron will be just enough to excite the second elec- 
tron. Hence, the momentum it will have subsequently 
to the collision will be that acquired from the field at 
the instant of rescattering. In terms of momentum-space 
constraints, this means that the region in the parallel 
momentum plane will collapse around piy = —A(t'). 

For a monochromatic field, A(t') — ±2^JU\ ) and 
E^(t',t") = 3.17C/ P HI- However, for a few-cycle 
pulse, —A(t') and E™^ x (t',t") will depend on the rescat- 
tering event within the pulse. This means that each spe- 
cific rescattering event will lead to a region in the Pi||P2|| 
plane, whose extension and center depend on the elec- 
tron return time and kinetic energy for a particular set 
of orbits. 



Similarly, Eq. ([9]) can be rewritten as 

{ Pn + A(t)} 2 = -2E { 2 e) - p 2 2± . (11) 

The above-stated expression shows that there will be a 
large drop in the yield with increasing transverse mo- 
mentum, as effectively widens the potential barrier 
through which the second electron must tunnel. The 
electron tunnels most probably at the laser field max- 
ima, for which A(t) <C 1 and p 2 \\ — 0. Hence, the 
momentum-space conditions for which ionization of the 
second electron is most probable lie approximately at 
(P2\\tP2x) = (0,0). For a few-cycle pulse, depending on 
the specific tunnel ionization event, the barrier will be 
narrower or wider. 

In summary, the constraints upon the momentum of 
the second electron determine the position of the cross- 
shaped distributions at the axes p n || = in the parallel- 
momentum plane, while those related to the momentum 
of the first electron determine the regions along these 
axes that will be populated, i.e., whether they will be 
long or short. Clearly, for a few-cycle pulse this length 
will vary with each particular event. 



III. QUANTUM-ORBIT ANALYSIS 

In this section, we will perform a quantum-orbit anal- 
ysis of the problem, with emphasis on how the dominant 
sets of orbits change with the CEP. The concept of "quan- 
tum orbits" is based on the fact that the solutions of the 
saddle-point equations can be related to the classical or- 
bits of an electron in a field and, still, retain information 
on quantum aspects such as tunneling and interference 
(for a broad overview see [lU). 

We employ the linearly polarized few-cycle pulse 
E(t) = —dtA(t)/dt, where the vector potential A(t) is 
given by 

A(t) = 2^ sin 2 (j^J sin(0 + ut)e,. (12) 

In Eq. (fT2")l , N denotes the number of cycles in the pulse, 
to is the field frequency, <f> the carrier-envelope phase, 
Up = E 2 /(Aoj 2 ) is the ponderomotive energy, E is the 
field amplitude and e z the polarization vector. Through- 
out, we choose the number of cycles as N = 4.3 and 
4> = 0i — 0O: where 0o = 60° is an offset value. This 
is well within the parameter range employed in experi- 
ments [l5|]- Throughout, we will refer to the phase (f>\ 
without the offset value to facilitate a comparison with 
the existing literature. The bound-state energies taken 
correspond to argon. Initially, both electrons are bound 
in the 3p state. The first electron recollides with the core, 
exciting a second electron to the 4s state. 
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FIG. 1: Schematic representation of the electric field E(t) 
and the corresponding vector potential A(t), for a short pulse 
of N = 4.3 cycles, whose shape is defined by Eq. ()12[) and 
whose carrier-envelope phases are the same as in the previous 
figure. The arrows indicate the approximate classical times 
around which the first electron leaves, in case it returns at a 
field crossing. The pairs of orbits are indicated by the labels 
Pair n, where n ranges from 1 to 8. The fields have been 
normalized to E(t)/Eo and A(t)/Ao, where Eq, Aq denote 
the field amplitudes. The carrier envelope phases are fa — 
65° = 5° ), fa = 155° (f/> = 95°), fa = 245° {(j> = 185°) and 
fa = 335° = 275°). 



A. Approximate ionization and rescattering times 

We will now identify the relevant sets of orbits for the 
pulse (|12l) . In Fig.[TJ we indicate the approximate ioniza- 
tion and rescattering times for the first electron, for the 
values of the CEP employed in this article. These sets 
of times are associated with the real part of the complex 
times t" and t', obtained from the solutions of the saddle- 
point equations ([f)])-©. Such solutions always occur in 
pairs, which, physically, correspond to the fact that the 
first electron may return along a shorter and a longer or- 
bit. In the figure, one may identify up to eight pairs of 
orbits, whose ionization and rescattering times vary with 
the CEP. In the following, we will refer to these pairs as 
Pair ro(ei), with n = 1, ...,8. For each pair, the electron 
will leave most probably at a local maximum and return 
most probably at the subsequent crossing. The most rel- 
evant pairs are those near the center of the pulse, as the 
field intensity in this case is higher. 

Fig. [2] shows the approximate times at which the sec- 
ond electron tunnels. These times are located around the 
field maxima, and also strongly depend on the CEP. The 
corresponding orbits will be referred to as Orbit nfa), 
with n ranging from 1 to 8. Once more, the orbits for 
which tunneling is expected to be most prominent are 
close to the center of the pulse. Note that the second elec- 
tron cannot tunnel from an excited state before the first 



FIG. 2: Schematic representation of the same electric field 
E(t) and the corresponding vector potential A(t) as in the 
previous figures, but highlighting the times around which ion- 
ization of the second electron is expected to occur, indicated 
by the squares in the figure. The orbits are indicated by the 
labels Orbit n, where n ranges from 1 to 8. The fields have 
been normalized to E(t)/Eo and A(t)/Ao, where Eq, Aq de- 
note the field amplitudes. 



electron rescatters. This implies that, for instance, ion- 
ization related to Orbit 2{ez) can only be caused by Pair 
l(ei), or ionization related to Orbit 3(e2) by Pairs l(ei) 
and 2(ei). Even though this sounds obvious, a rigorous 
treatment of causality can be a non-trivial issue, espe- 
cially within the context of the steepest descent method, 
and requires an extensive modification of the contours to 
be taken into account. For a detailed discussion, see our 
previous work [27]. 



B. Solutions of the saddle-point equations 

1. First electron 

We will now bring together the intuitive picture dis- 
cussed above and the solutions of the saddle-point equa- 
tions. We will start by focusing on the first electron, for 
which the saddle-point equations ([B])-© give the com- 
plex ionization and rescattering times. We will study 
Pairs 3(ei), 4(ei ) and 5(ei), as, for a wide range of CEPs, 
such pairs are expected to lead to the most relevant con- 
tributions for the specific pulse chosen. For the cases 
studied in this section, the saddle-point solutions have 
been obtained for vanishing momenta p\±^ and different 
CEP values. 

Fig. [3] displays the real and imaginary parts of the 
above-mentioned solutions, for Pair 3(ei). As an over- 
all feature, the real parts of the ionization and return 
times are centered around = —A(t'), which is in the 
negative parallel momentum region (see Figs G2a) and 
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FIG. 3: Real and imaginary parts of the complex ionization 
and rescattering times t" and t' of the first electron obtained 
from the solutions of the saddle-point Equations ©-([H]) cor- 
responding to Pair 3 in [T] for vanishing transverse momen- 
tum pi± — 0. Panels (a) and (b) give the real parts of such 
solutions, and panels (c) and (d) the imaginary parts. The 
bound-state energies correspond to argon assuming a 3p — > 4s 
excitation (Ei g — 0.58 a.u., Z?2 fl = 1-02 a.u. and E2 e = 0.40 
a.u.) in a few-cycle pulse of frequency u) = 0.06085 a.u. and 
intensity I = 3 x 10 14 W/cm 2 . 



(b), respectively). This is the most probable momen- 
tum for the first electron upon rescattering. Further- 
more, Re[wi"] and Refwt'] almost coalesce at two spe- 
cific values of the parallel momentum components pin. 
These are the minimum and the maximum momentum 
values for which rescattering exhibits a classical counter- 
part. Outside this region, this process is forbidden and 
the corresponding transition amplitude is exponentially 
decaying. 

The figure shows that the ionization and rescatter- 
ing times, the center and the extension of the classi- 
cally allowed region strongly depend on the CEP. In 
particular, there is a very good agreement between the 
center of the region shown in Figs. EJa) and (b) and 
the ionization and return times indicated by arrows in 
Fig. [TJ which were identified employing classical argu- 
ments. As the CEP increases, the real parts of such 
times move towards the pulse turn-on. This will lead 
to a decrease in the extension of the classically allowed 
region from —\^JTT P < pu < Q.7^/Up for fa = 65° to 
-2^/TT v < p!|| < -0.8 v /0p for fa = 335°. This can 
be physically understood by direct inspection of Fig. [IJ 
and by bearing in mind that the classical limit of the 
kinetic energy of the first electron upon return is given 
by Ef& = (A(t') - A(t")) 2 /2^. For fa = 65°, at the 
ionization time related to Pair 3(ei), A(t") ~ 0, and, at 
the rescattering time, the vector potential A(t') is around 



its peak value Ao. Hence, the kinetic energy of the first 
electron upon return is high. This implies that, for this 
phase, the contributions of Pair 3(ei) will populate a 
large momentum region. As the CEP is increased up 
to fa — 335°, the instantaneous vector potential A(t') 
decreases to less than O.8A0. Furthermore, due to the 
lack of monochromaticity of the field, A(t") is no longer- 
vanishing near the pulse turn-on. Both effects lead to 
a decrease in this momentum region. Apart from that, 
the center of the classically allowed region, which is de- 
termined by the most probable momentum the electron 
may have upon return, moves from roughly — 2^/Up to 
— lAy/Up. This is due to the fact that, close to the peak 
of the pulse, the former value is a good approximation 
for — A(t'), while, near the edges of the pulse, its lack of 
monochromaticity plays an increasingly important role. 

The imaginary parts of t" and t', shown in the remain- 
ing panels of Fig. |3l give valuable information on whether 
the process in question is allowed or forbidden, or on 
the overall probability related to a specific process, e.g., 
tunneling ionization. For instance, Imfcjf"], displayed in 
Fig. EJd), sheds some light on how the width of the po- 
tential barrier that the electron must overcome in order 
to reach the continuum varies for Pair 3(ei), with regard 
to the CEP. The larger Imfwt"] is, the wider is the po- 
tential barrier through which the electron must tunnel 
[34| . The picture shows a marked increase in Im[i"] as 
the CEP varies from fa = 65° to fa = 335°. This can 
be understood with the help of Fig. Q] According to this 
figure, for fa = 65°, the first electron tunnels near a lo- 
cal maximum for which the instantaneous electric field is 
O.SEq. As the CEP increases, the local maximum associ- 
ated with tunnel ionization for Pair 3(ei) decreases down 
to less than 0.4-Eo for fa = 335°. Hence, the potential 
barrier widens and the contributions of this specific pair 
become less and less relevant. Note that Im[u;i"] ^ 
throughout, as tunneling is not classically allowed. 

Finally, Fig. [3jd) shows the behavior of Im[wf'], which 
is associated to the rescattering time, with regard to the 
CEP. In contrast to what happens to Im[wt"], this imag- 
inary part vanishes at the momentum range between the 
values of pu for which Re[wt"] and Re[wi'] almost co- 
alesce. This is related to the fact that, in this region, 
rescattering is classically allowed. This region decreases 
in extension for increasing values of fa . Once more, this 
reflects the fact that Pair 3(ei) loses relevance. 

The above-mentioned loss of relevance does not happen 
to all pairs of orbits, but will depend strongly on how the 
corresponding start and return times are located within 
the pulse. In fact, it may occur that, as the CEP in- 
creases, a specific set of trajectories becomes dominant, 
or even remains relatively stable. In order to understand 
this issue, it suffices to analyze the real and imaginary 
parts of the ionization time t" . An example is provided 
in Fig. |4]for Pair 4(ei) and Pair 5(ei). The contributions 
of the former pair to the NSDI distributions remain rel- 
atively stable, while the latter pair increases in relevance 
within the CEP range studied. The real and imaginary 
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FIG. 4: Real and imaginary parts of the complex ionization 
times t" of the first electron obtained from the solutions of the 
saddle-point Equations ©-([8]) corresponding to Pairs 4 and 5 
in [U for vanishing transverse momentum pj± — 0. Panels (a) 
and (b) give the real parts of such solutions, and panels (c) 
and (d) the imaginary parts. The field and atomic parameters 
are the same as in the previous figures. 
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FIG. 5: Real and imaginary parts of the complex ionization 
times t of the second electron obtained from the solutions of 
the saddle-point Equation ([9]l) corresponding to Orbits 4 to 
7 in [2] for vanishing transverse momentum p2± = 0. Panels 
(a) to (d) give the real parts of such solutions, while panels 
(e) to (h) depict the imaginary parts. 



Pair 5(ei). Fig.|4]Jd) shows that such imaginary parts are 
comparable to those observed for Pair 4(ei) for a wide 
range of phases. Hence, we expect both sets of orbits to 
compete, until Pair 5(ei) becomes dominant. 



parts of the complex ionization time t" are depicted in 
Fig. S] for both orbits. 

Fig. SJa) shows that, for Pair 4(ei), the classically al- 
lowed region, centered at a positive parallel momentum, 
remains roughly the same throughout. This happens be- 
cause the absolute value of the vector potential A(t'), 
and hence the kinetic energy obtained by the electron 
upon return, remains quite large for the CEP range stud- 
ied. The behavior of Pair 5(ei), however, displayed in 
Fig.0Jb), is markedly different. For this pair, the exten- 
sion of the classically allowed region increases substan- 
tially as the CEP increases. This is expected as, in this 
case, the approximate ionization and rescattering times 
in Fig. [1] move from the turn off to the center of the pulse 
with increasing CEP. Consequently, the vector potential 
A{t') at the instant of rescattering increases from approx- 
imately 0.4Ao to Aq. This leads to a substantial increase 
in the electron kinetic energy upon return. 

Another important issue determining the dominance 
of a pair of orbits is the width of the potential barrier 
through which the first electron tunnels, which can be 
roughly inferred from Im[aj£"]. For Pair 4(ei), Im[uit"} re- 
mains relatively stable throughout, as shown in Fig. HKc) . 
An inspection of Fig. Q] shows that, indeed, the instanta- 
neous electric field strength \E(t")\ related to Pair 4(ei) 
is rather large. Hence, one expects the corresponding po- 
tential barrier to be quite narrow. For increasing CEP, 
Im[i"] increases slightly. This is due to the fact that 
Ke[ut"} related to this pair moves towards the pulse turn 
on and the potential barrier at the ionization time t" be- 
comes slightly wider. A similar analysis can be performed 
in the imaginary parts of the saddle-point solutions for 



2. Second electron 

For the second electron, a set of orbits will only be 
dominant if it is excited in such a way that its subsequent 
tunneling probability is high. In order to obtain this in- 
formation, we will look at the tunnel ionization times of 
the second electron, which can be obtained by solving 
the saddle-point equation ([9]). The real and imaginary 
parts of these solutions are displayed in Fig. [5] (upper 
and lower panels, respectively), for vanishing transverse 
momentum p 2 ±_ = 0. We only consider Orbits n(e 2 ) for 
n ranging from 4 to 7. The preceding orbits, i.e., Or- 
bits n(e 2 ) with n = 1,2,3, will not play an important 
role, as they violate causality with regard to Pair 3(ei). 
Note that, for the second electron, the orbits are well 
separated for all momentum ranges. There is also no 
classically allowed region, as tunneling is an intrinsically 
quantum mechanical phenomenon. Hence, Im[oji] ^ 
throughout. Once more, one may draw an intuitive pic- 
ture by relating the solutions in Fig. [5] to the simplified 
arguments illustrated in Fig. [2] 

For all phases, we observe that, for parallel momen- 
tum p 2 \\ — 0, Refwi] corresponds to the times for which 
the pulse has local cxtrema, marked by the rectangles in 
Fig.[2j These extrema vary from the center until the pulse 
turn off, and correspond to the most probable ionization 
times for the second electron, i.e., when the potential bar- 
rier is narrowest for a specific orbit. A maximal ioniza- 
tion probability near p 2 \\ =0 and its relation to the nar- 
rowest potential barrier can also be inferred by inspecting 
Im[a;i], which exhibits a minimum around this value (see 
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lower panels in the figure). The above-stated features, 
i.e., the electron tunneling most probably at peak fields 
and a corresponding minimum at Im[t] at such times, are 
also observed for a monochromatic field (23[. 

However, due to the lack of monochromaticity in the 
pulse, there are sometimes small deviations from p 2 || = 0, 
especially if the orbits are near the pulse turn on and 
off. This feature can be seen very clearly, for instance, in 
Fig. EJh), which corresponds to a tunnel ionization event 
close to the pulse turn off, i.e., to Orbit 7(e2). Further- 
more, Im[u;i] is not symmetric with regard to its value 
at P2\\ — 0. This would be the case for a monochro- 
matic field, as all cycles, and hence the regions around 
each local field maximum, would be identical for a con- 
tinuous wave (23|. In a few-cycle pulse, however, one 
expects the potential barrier through which the second 
electron must tunnel to narrow, or to widen, as the ion- 
ization times approach or distance themselves from the 
center of the pulse, respectively. This leads to asymme- 
tries in Im[o;i] for positive and negative momenta. In 
the specific pulse discussed in this work, for Orbits 5(e2) 
and 7(e 2 ) [Figs. [5jf) and (h), respectively], Im[wi] de- 
creases for p 2 \\ > 0, in comparison to the region for which 
p 2 \\ < 0, while, for Orbit 6(e2), the opposite behavior 
holds. This asymmetry can be understood by inspect- 
ing the respective upper panels, together with Fig. [2] 
For Orbits 5(e2) and 7(e2), if the electron leaves before 
the corresponding field maxima, its momentum will be 
positive (see Figs. [Sfb) and (d), respectively). An in- 
spection of Fig. [2] shows that these temporal regions are 
located closer to the center of the pulse, compared to the 
time regions subsequent to the local field maxima. In 
contrast, for Orbit 6(e2), earlier ionization times, and a 
larger proximity to the central region of the pulse, corre- 
spond to negative momenta. For orbit 4(e2), the above- 
mentioned asymmetry varies, as shown in Fig.[5Je). This 
happens because, throughout, the ionization times are 
very close to the pulse center. Hence, whether ioniza- 
tion times prior or subsequent to the local maxima will 
correspond to higher or lower local field intensities will 
depend on the carrier-envelope phase. 

The dominant sets of orbits in a few cycle pulse can be 
identified by analyzing Fig. [5j As a general feature, as 
4>i increases, the real parts of the ionization times move 
towards lower values, as shown in Figs. [5](a)-(d). This 
is consistent with Fig. [5] and with the previous analy- 
sis, performed for the first electron. Some orbits, such 
as Orbits 5(e2) - 7(e2), move towards the center of the 
pulse, while others, such as Orbit 4(e2), move away from 
it. This behavior, and its consequences for the shape 
of the electron-momentum distributions, can be inferred 
by analyzing the imaginary parts of the ionization times, 
displayed in Figs. [S](e)-(h). These panels show that the 
most important orbits for the second electron will be Or- 
bits 4(e2) and 5(e2). This is due to the fact that, for 
these orbits, near p 2 \\ — 0, 0.5 < Im[a;i] < 0.75, while, 
for the remaining orbits, Im[u;i] is in general larger. This 
implies that the potential-energy barriers through which 



the electron must tunnel are narrower for Orbits 4(e2) 
and 5(e2). Interestingly, Im[u;i] remains stable for Or- 
bit 4(e2) throughout. This may be understood by in- 
specting the local maximum of E(t) in Fig. [5] related to 
this orbit. The figure shows that the instantaneous field 
strength at this maximum remains located near 0.8i?o for 
the CEP range considered. In contrast, for Orbit 5(e2), 
Imjwt] decreases systematically from around 0.75 to 0.5 
as the phase varies from from <fii = 65° to (pi = 335° (see 
Fig- EI 1 ))- This occurs because the instantaneous field 
strength \E(t)\ related to this orbit increases from less 
than 0.5i?o to almost the full amplitude Eq in this phase 
interval, as shown in Fig. [2] 

For the remaining pairs of orbits, there is always a de- 
crease in Im[a;t] with regard to increasing CEP values. 
This is due to the fact that the real parts of such times 
move from the pulse turn off towards the center of the 
pulse, so that the effective potential barrier becomes nar- 
rower. In particular, we expect vanishingly small contri- 
butions throughout for Orbit 7(e2), as in this case Im[u;i] 
is much larger than for the other orbits. 

IV. ELECTRON-MOMENTUM 
DISTRIBUTIONS 

In this section, we will apply the information gained 
from the quantum-orbit analysis in order to determine 
the regions in momentum space that should be occupied 
by the electron-momentum distributions. These predic- 
tions will be compared with the outcome of the actual 
computations, performed by setting Vp ie ^ g — const., 
Vp 2 e = const. This guarantees that both the excitation 
and ionization prefactors do not introduce any momen- 
tum bias, and the distributions represent the momentum 
constraints discussed in the previous section. 

A. Partial distributions 

For simplicity, we will commence by analyzing the par- 
tial momentum distributions 

F^( Pn )= J \M^( Pl )\ 2 d 2 Plx (13) 

and 

^ (2) b 2 ||) = J \M^(p 2 )\ 2 d 2 p 1± , (14) 

where the partial transition amplitudes for the first and 
the second electron read 

/oo /*oo /• 
dt" dt' d 3 kV £>l V (9) e iSl(pi,Kt",t') 
-oo Jt" J 

(15) 

and 

/oo 
dtV^e^^, (16) 
-OO 
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respectively and the transverse momentum components 
are integrated over. Note that, due to causality, the to- 
tal ionization probabilities are not the product of such 
functions 



27]. 



In Fig. [6l we display F^'(pi») and F^ 2 \p 2 \\) computed 
for each of the contributions described in the previous 
section. For the first electron, the figure illustrates very 
clearly the loss of relevance in Pair 3(ei) discussed above 
as the CEP increases (see Figs-UJa) and (b)). The contri- 
butions from such a pair, mostly located in the negative 
parallel momentum region, become vanishingly small al- 
ready for <t>\ = 245°. In fact, even for fa = 65° this pair 
is not dominant, as the partial probability density asso- 
ciated with Pair 4(ei), mostly in the positive momentum 
region, is higher. This is a consequence of the fact that 
the field strength at the corresponding ionization time is 
higher for Pair 4(ei), so that the first electron tunnels 
through a broader barrier (see Fig. [T]). In Figs. Eth) to 
(d), the partial momentum distributions also show that 
Pair 5(ei) is very important over a large CEP range, and 
in fact provides the dominant ionization pathway for the 
first electron for fa = 155° and fa = 245° (see Figs, fflb) 
and (c)). This is in agreement with with Fig. [TJ which 
shows that, for a large range of phases, the instantaneous 
field at ionization is near its absolute maximum for this 
pair. In fact, only for fa = 345° = 335° does Pair 6(ei) 
become dominant. Another interesting feature is that 
the curves corresponding to a specific Pair n(ei) are the 
mirror images of those related to Pair n + l(ei) if the 
CEP is shifted in A0 = 180°. This holds for all partial 
probabilities displayed in Figs.|6][a) to (d). For instance, 
for Pair3(ei) in Fig.^a) is equal to F^(-p n ) 
for Pair 4(ei) in Fig. [Bfc) and so on. Finally, the peaks 
of the partial distributions agree with the predictions in 
Sec. UTT] obtained from the solutions of the saddle-point 
equations. 



For the second electron, the partial probabilities 
F( 2 )(p 2 ||) confirm that the main ionization channel is 
via Orbit 4(ei) for a large range of CEPs, as shown in 
Figs. [Hie) and (f). The contributions from Orbit 5(e2) 
only start to compete with this channel at fa — 245°, 
and eventually become dominant at fa = 335° [Figs.EJg) 
and (h), respectively]. One should note that the symme- 
try F^ 2 \p 2 \\) for Orbit n(e 2 ) corresponds to F^ 2 >{— f> 2 ||) 
for Orbit n+ l(e2) also occurs upon a phase shift Afa = 
180°. In the figure, however, it cannot be observed as Or- 
bit 3(e2) is missing. This orbit corresponds to ionization 
events triggered by Pairs l(ei) or 2(ei), as illustrated 
in Figs. [T] and [5J Such pairs are irrelevant due to the 
weak field amplitudes involved [3^]. In Table U we pro- 
vide a summary of the most relevant orbits encountered 
for the first and second electron, in decreasing order of 
importance. 
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FIG. 6: Partial electron momentum distributions F^(p n \\) 
(n — 1, 2) fEqs,U3land ll4[) computed using the three most rel- 
evant individual ionization pathways for the first and second 
electron according to Figs. Q] and [5] and Table U for the same 
field and atomic parameters as in the previous figure. The up- 
per and lower panels refer to the first and second electron, re- 
spectively. From left to right, the phases fa = 65° [panels (a) 
and (d)], fa = 155° [panels (b) and (e)], fa = 245° [panels (c) 
and (f)] and fa = 335° [panels (d) and (g)] have been taken. 
The probability densities have been normalized to slightly be- 
low unity for the dominant contributions in each panel, in or- 
der to facilitate a more direct comparison. The dotted lines 
in panels (a) to (d) indicate the momenta i:2s/Up, where the 
partial distributions related to a monochromatic driving field 
are expected to be peaked. 



B. Correlated electron momentum distributions 

In view of the discussion presented above, one could 
conclude that a pair of orbits will lead to dominant con- 
tributions in the electron-momentum distributions if (i) 
the corresponding potential barrier through which the 
electrons tunnel is as narrow as possible, as this will lead 
to a high ionization probability; (ii) the kinetic energy of 
the first electron upon return is as high as possible, as 
this will populate a large region in momentum space. One 
should bear in mind, however, that RESI is a correlated 
two-electron process. This implies that the shapes and 
regions populated by the NSDI electron momentum dis- 
tributions will be determined by the interplay between 
the dominant contributions from the first and second 
electron. 

In order to understand this issue, we compute the cor- 
related electron- momentum distributions 

F{pi\\,p 2 \\) = J J d 2 p 1± d 2 p 2 ±\M(p 1 ,p 2 ) + p 1 op 2 | 2 , 

(17) 

as functions of the momentum components p n » (n = 1,2) 
parallel to the laser- field polarization. In Eq. (fTTj) . un- 
less otherwise stated, M(pi,p2) is the transition ampli- 
tude (P) associated with the coherent sum of all processes 
over all sets of orbits, both for the first and the second 
electron. These transition amplitudes are symmetrized 
to account for the fact that both electrons are indistin- 
guishable. The transverse-momentum components are 
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(f>i (degrees) 


First electron Second electron 


65° 


1. Pair 4 (+) 1- Orbit 4 

2. Pair 3 (-) 2. Orbit 5 
Pair 5 (-) 3. Orbit 6 


155° 


1. Pair 5 (-) 1. Orbit 4 

2. Pair 4 (+) 2. Orbit 5 

3. Pair 3 (-) 3. Orbit 6 


245° 


1. Pair 5 (-) 1. Orbit 4 

2. Pair 4 (+) 0rbit 5 
Pair 6 (+) 2. Orbit 6 


335° 


1. Pair 6 (+) 1. Orbit 5 

2. Pair 5 (-) 2. Orbit 4 

3. Pair 4 (+) 3. Orbit 6 



TABLE I: Most relevant orbits for the first and second elec- 
trons for the values of the carrier-envelope phase employed in 
this work, in order of decreasing importance. A single num- 
ber for more than one orbit indicates that their contributions 
are comparable or competing, while different numbers indi- 
cate the clear dominance of a pair for the first electron, or an 
orbit for the second electron. The signs (±) indicate whether 
a specific pair of orbits leads to electron momentum distri- 
butions I-F^Hpih)! 2 peaked at positive (+) or negative ( — ) 
momentum. Subsequently, this will imply that the contribu- 
tions triggered by such a pair of orbits will populate either he 
positive or the negative half axis in the parallel-momentum 
plane. 



integrated over. These distributions are plotted in Fig. [7] 

For <pi = 65°, the distributions occupies a broad region 
along the negative half axes = (jpju < 0), with i ^ j, 
as displayed in Fig. [TJa) . This illustrates the importance 
of Pair 3(ei), as a direct comparison with Fig.[7][e) shows. 
Even if this pair is not related to the most prominent ion- 
ization event for the second electron, upon rescattering, 
it triggers the dominant tunnel ionization channel for the 
second electron, i.e., along Orbit 4(e2). This counterbal- 
ances the influence of Pair 4(ei), which, according to the 
partial distributions in Fig. [5]and our previous line of ar- 
gument, is the prevalent ionization pathway for the first 
electron for this specific CEP value. Contributions from 
Pair 4(ei), along the parallel momenta positive half axes 
are also present. These contributions are, however, com- 
parable, or even slightly weaker, as rescattering along 
Pair 4(ei) can only lead to ionization along Orbit 5(e2), 
Orbit 6(e2) or even later orbits. Furthermore, because 
an electron returning along Pair 3(ei) acquires a higher 
kinetic energy than if it returns along Pair 4(ei), the 
momentum region populated by events related to former 
pair is larger. This fact is also observed for the par- 
tial momentum distribution in Fig. H]Ja), which exhibits 
comparable probability densities over a broad momentum 
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FIG. 7: Correlated RESI electron momentum distributions 
as functions of the parallel momentum components pjri, p 2 || 
for the same atomic and field parameters as in the previous 
figures. In panels (a) to (1), constant ionization and excita- 
tion prefactors V pie ,kg, V P2e prefactors have been employed, 
while in panels (m) to (p) the prefactors related to the exci- 
tation from a 3p to a 4s state by a contact type interaction 
have been used. The explicit expressions for these prefactors 
is given in [2l|. From left to right, the carrier-envelope phases 
are <j>i = 65° [panels (a), (e), (i) and (m)], <j>\ = 155° [pan- 
els (b), (f), (j) and (n)], 0i = 245° [panels (c), (g), (k) and 
(o)], 4>i = 335° [panels (d), (h), (1) and (p)]. When constant 
prefactors are used, from the top to the bottom of the figure, 
we plot the distributions obtained as follows. Panels (a) to 
(d): all pairs of orbits for the first electron, and all possible 
ionization pathways for the second electron displayed in Ta- 
ble U Panels (e) to (h): all ionization channels for the second 
electron specified in Table U and Pair 3(ei) [panel (e)], 4(el) 
[panels (f) and (g)] and 5(ei) [panel (h)] for the first electron. 
Panels (i) to (1): all pairs of orbits for the first electron speci- 
fied in Table[T] but only the ionization pathways for the second 
electron immediately after that specific ionization event, i.e., 
Orbit 4(e 2 ) for Pair 3(ei), Orbit 5(e 2 ) for Pair 4(ei) and Orbit 
6(ei) for Pair 5(ei). In panels (m) to (o), we have included 
all orbits for the first and second electron specified in Tabled] 
For panels (a) to (1), the color scales range from to 1 x 10 , 
while in panels (m) to (o) they range from to 1 x 10~ 20 . The 
white lines in the figure indicate the antidiagonal pu\ — — Pall • 



region. 

For <p = 155°, in contrast, Pair 4(ei) prevails (see 
Fig. [TJb) as compared to [T^f)). This leads to a shift 
in the distribution towards the positive half axis of the 
Pi||P2|| plane. Physically, this can be attributed to the 
loss of relevance related to Pair 3(ei), together with the 
fact that the ionization channel along Orbit 5(e2) be- 
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comes more prominent. Hence, ionization along Orbit 
4(e2) can no longer counterweight the other effects. This 
is in agreement with the previous discussions in Sec. ITTT1 
and Fig. [5] Note, however, that the partial probability 
density F^{pu\) associated with Pair 5(ei) is dominant 
in Fig. EJb). Nevertheless, Pair 5(ei) can only lead to 
ionization events related to Orbit 6(e2) or later. These 
events are too close to the pulse turn off to play a signif- 
icant role. 

The distributions obtained for <\>\ — 245° [Fig. [7Jc)], 
on the other hand, show that, as the CEP increases, Pair 
5(ei) becomes more relevant and starts to influence the 
overall distributions. For this specific phase, the ioniza- 
tion channel for first electron along Pair 3(ei) is neg- 
ligible, regardless of the subsequent ionization events. 
Ionization along Pair 4(ei) is relatively small. It may, 
however, cause ionization of the second electron along 
Orbit 5(ei), which is quite prominent. Pair 5(ei) is the 
most prominent ionization channel for the first electron, 
but may lead to ionization only along Orbit 6(e2) or at 
later times. Hence, Pair 4(ei), whose contributions are 
presented in Fig. UJg), still determines the momentum 
regions to be occupied. Interestingly, this distribution 
is the mirror image of that obtained for (f>i = 65° with 
regard to (p 1 \\,p 2 \\) -> (~Pi\\, ~P2\\)- 

Finally, for (f>i — 335°, displayed in Fig. [TJd), the dis- 
tributions are almost entirely concentrated along the neg- 
ative half axis p^\ = (pj» < 0), with i ^ j. This is a 
consequence of the fact that, together, the pathways re- 
lated to Pair 5(ei) for the first electron, and Orbit 6(e2) 
for the second electron determine the momentum regions 
to be populated (see Fig. [7th) for comparison). A similar 
interpretation provided when discussing the distributions 
obtained in Fig. [7ta), for <fii = 155°, applies, with the 
difference that all the indices must be shifted by 2, i.e., 
instead of Pairs 3(ei) and 4(ei), now one must consider 
Pairs 5(ei) and 6(ei) and the subsequent orbits for the 
second electron. Note once more that this distribution is 
the mirror image of that depicted in Fig. UJb). 

Apart from the effects discussed above, throughout, 
there exist residual fringes parallel to the axis p n ii = 0. 
These fringes are related to the fact that, in our com- 
putations, we have included all the relevant orbits along 
which the second electron may tunnel, subsequently to 
being excited by a particular pair of orbits. Quantum 
mechanically, the transition amplitudes related to these 
pathways interfere and even partly survive the integra- 
tion over the transverse momentum components. In a 
more realistic scenario, however, we expect these fringes 
to be absent, due to the fact that the excited bound state 
from which the second electron is released is strongly de- 
pleted. Hence, ionization would mainly occur around the 
field maximum closest to the time of excitation. 

In order to mimic depletion, in Figs.[7Ji) to (1) we con- 
sider only this ionization pathway. As an overall feature, 
the above-mentioned structure disappears. Furthermore, 
the correlated probability densities are now slightly dis- 
placed from the axes p n \\ = 0. This occurs because, 



due the lack of monochromaticity of the few-cycle pulse, 
Im[wi] no longer exhibits a minimum at p^ =0 (see dis- 
cussion in Sec. lIII[) . If many ionization events are consid- 
ered, this effect tends to average out and this asymmetry 
is less prominent. Concretely, the upper parts of the dis- 
tributions in Fig.[7Jare slightly displaced towards the first 
quadrant of the parallel momentum place. This is caused 
by the fact that ionization along Orbit 5(e2), which is the 
main ionization channel for the second electron related 
to Pairs 4(ei) and 6(ei), is favored for positive momenta 
(see Fig. |UJf)) and subsequent discussion). In contrast, 
the lower parts of such distributions are slightly shifted 
towards the third quadrant. This is due to the fact that 
the main ionization channel for the second electron in this 
momentum region is either along Orbit 4(e2) (4>i = 65° 
and fa = 155°) or along Orbit 6(e 2 ) (fa = 245° and 
fa = 335°). For both orbits, ionization is slightly en- 
hanced for negative momenta, as shown in Figs.Sfe) and 
(g). Finally, in Figs.[7Jm) to (p) we incorporate the pref- 
actors Vp l9l kg and V P2e corresponding to an excitation 
process from the 3p to the 4s state a contact-type inter- 
action 1^2 (i"i, £2) = S(ri — r 2 ). For the explicit expres- 
sions see Ref. [21(. In general, these prefactors introduce 
a bias towards low momenta, so that the distributions are 
much more focused around the origin p^ = p 2 \\ = 0. For 
a Coulomb- type interaction V12 = l/|i"2 — we have 
verified that even lower momenta are favored. 



V. CONCLUSIONS 

In this work, we highlight the influence of the carrier- 
envelope phase (CEP) for the recollision-excitation with 
subsequent tunneling ionization (RESI) pathway in NSDI 
with few-cycle pulses. The electron-momentum distribu- 
tions are quite sensitive with regard to the CEP, and, as 
this parameter varies, move from the region below to the 
region above the anti-diagonal pi\\ = —p 2 \\ i or y i ce versa, 
in the plane spanned by the electron-momentum compo- 
nents parallel to the laser field polarization. Similarly to 
our previous studies in 0, H| performed for the simpler 
NSDI mechanism of electron-impact ionization, we show 
that all features encountered can be explained in terms of 
electron trajectories. Indeed, a detailed analysis of where 
these trajectories are located within the pulse, and, in the 
context of the steepest descent method, of the real and 
imaginary parts of the ionization times of both electrons 
and the rescattering time of the first electron provide 
a consistent picture related to the shapes, maxima and 
regions populated by the RESI electron-momentum dis- 
tributions. All such features are determined by the in- 
terplay between the dominant sets of trajectories for the 
first and second electron. 

Qualitatively, the results in the present publication 
agree with those observed experimentally in [15| , in which 
it has been reported that, depending on the CEP, the 
electron-momentum distributions shift from the region 
below to the region above the anti-diagonal pi|| = ~P2\\- 
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Quantitatively, however, we find a much larger probabil- 
ity density around the origin pjii = p 2 \\ = than those 
encountered in [l5| ■ This is particularly true if we con- 
sider the prefactors V P2e , V Vie ^ g associated with an ex- 
citation from the 3p to the 4s state in Argon. In this lat- 
ter case, the distributions become much more localized 
around the origin. These discrepancies may be related 
to several issues. First, there is always an uncertainty 
in the experimentally measured peak intensities. If this 
intensity is smaller than that employed in our computa- 
tions, this means that the region around = p 2 \\ = 
is less populated than in Fig. [JJ (see [2l|, |23| | for details). 
Second, in a more realistic scenario, it may be that other 
excitation processes than that considered in this work 
play an important role. This implies that the features 
specific to particular bound states may average out, and 
one may approach the results obtained when the pref- 
actors were left out. Third, the momentum constraints 
related to the first and the second electron, which have 
been derived within the context of the strong-field ap- 
proximation [2l|, (27j neglect the presence of the binding 
potential when the electron is in the continuum, and as 
well bound-state depletion. Both issues may affect such 
constraints. 



Finally, we would like to comment on the role of deple- 
tion on the RESI electron-momentum distributions. In 
[HI ], the probability densities at nonvanishing momenta 
which led to the axes of the cross have been attributed 
to the second electron. It has been argued that, due to 
depletion, the second electron left before the peak field 
with non-vanishing momentum. Our results, together 
with the constraints stated in this paper, suggest that 
depletion will mainly shift the probability density with 
regard to the axis p„ii = (see discussion of panels (i) 
to (1) in Fig. [7]). This effect, however, seems to be rel- 
atively small. According to our model, the shift in the 
probability density away from the origin pi\\ = p 2 \\ is 
mainly caused by the first electron, which, upon recolli- 
sion, acquires the additional momentum — A(t) from the 
field according to the constraints stated in Sec. [Til The 
issue of depletion is not yet fully understood, and will be 
pursued in future work. 
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